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We study the dynamics of the transmission of the hanta virus infection among mouse populations, 
taking into account, simultaneously, seasonal variations of the environment and interactions within 
two classes in the mouse population: adults and subadults. The interactions considered are not 
symmetric between the two age-organized classes and are responsible for driving the younger 
members away from home ranges. We consider the case of a bounded habitat affected by seasonal 
variations. 



PACS numbers: 87.10.Ed, 87.23. Cc 



INTRODUCTION 

Theoretical investigations of population dynamics of 
animals such as mice or other rodents derive their current 
importance both from their direct relevance to the spread 
of epidemics such as the Hantavirus [l|, 0, 0, 3 [H| but also 
from the interplay of the approaches of physics and ecol- 
ogy that such investigations encourage 
Our interest in the present paper is on prevalent ef- 
fects of the competitive interaction among adults and 
subadults among rodent populations. Several ecological 
studies have shown the importance of such an interaction 
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l3|. The interaction can be studied by considering 
an age-structured population of mice, i.e., by splitting 
the population into classes organized by age. There is 
competitive struggle for territory associated with the fact 
that the occurrence of infected cases of several types of 
Hantavirus is mostly observed among adult individuals. 
This, in turn, is related to the proposed mechanism of 
infection among rodents: territorial fights that produce 
wounds on the animals as a tangible and verifiable conse- 
quence. The subadult population, consisting of smaller- 
sized individuals, tends to maintain its distance from the 
fights, and is thus less susceptible to infection. The effect 
of remaining not too close to the source of the fights, viz., 
the adult population, subadult individuals are forced to 
abandon the already colonized spaces and are therefore 
driven to new habitats, most of the time less suitable for 
survival. 

We propose here a model considering an age-structured 
population composed by the two well differentiated 
groups we have mentioned above: subadults and adults. 
A third group, composed by juvenile individuals may be 
taken into account but will be left out of our calculations 
as relatively unimportant. 

Like previous models 0, 0, 0, S], ours is based on a set 



of Fisher-like equations [13| to describe the evolution of 
the population of individuals of each group whose densi- 
ties we will respectively denote by Ma and My. Here, the 
suffix a stands for adults and the suffix y for the young 
individuals, the subadults. 



THE MODEL 

The specific set of equations in our model would be, in 
the absence of territorial interactions, 
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Here we are considering different diffusion coefficients Da 
and Dy for the adults and the subadults, to take into 
account the fact that adult individuals tend to remain 
in a rather bounded area, considered as each individual 
territory, while the subadults tend to move further and 
probably faster. One would therefore expect Da < Dy. 
The term in each of the above equations which is associ- 
ated with the competition for resources, is characterized 
by what is called the environmental parameter K and is 
proportional to the so-called carrying capacity. We will 
take K to be time dependent because of seasonal vari- 
ations. We denote the total population of the mice by 
M = Ma + My. The rate /i is one of transformation 
of subadults into adults through the process of matu- 
rity, thus is proportional to the subadult population and 
appears with reversed sign in the two equations. The 
subadults are born from the adults at (3, the birth rate. 
Death of the adults occurs at rate S and we have omitted 
it from the subadult population equation for simplicity 
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and to reflect the simplifying assumption that natural 
death visits only the adults. If predators were introduced 
into our considerations different rates might be put into 
the equations including a death rate for the subadults. 

In this work we propose to modify the above set of 
equations to take into account the ecologically impor- 
tant experimental observation related to territorial flghts. 
Due to competition for the conquest and preservation of 
the home range, adult mice tend to fight among them- 
selves. Indeed, these fights have been suspected to con- 
stitute one of the most important ways of transmission 
of the Hantavirus among rodents of the same species. 
As a result of these territorial threats of the adults, the 
subadults, being smaller in size, are forced to abandon 
already colonized environments and move towards unoc- 
cupied spaces. There are at least two different ways to 
include this tendency into the equation. To take into 
account the fact that subadults will tend to move away 
from places with high adult occupation, we may assume 
the subadult fiux to be proportional to the gradient of the 
adult population, VM^, at the same time, no subadult 
flux is possible in the absence of young individuals. Thus, 
the interaction term should be proportional to the pop- 
ulation of the subadults My as well as to the adult gra- 
dient. The resulting set of equations is 



dMg 
dt 



dMy 
dt 
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with K a constant that determines the strength of the 
interaction. 

The other manner that young individuals might tend 
to get away from populated areas is by creating a flux of 
their own class which is proportional to the local density 
of the other class, i.e., the adult. Under these conditions 
we will have 
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dt 
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Comparison of the two new interaction expressions. 



and 



V{MyVMa) = {VMa){VMy) + MyVH'Ia 



V{MaVMy) = {VMa)(VMy) + MaVH'Iy 



shows that they share a common term proportional to the 
gradient of each of the classes and another term which 



in one case is proportional to the density of subadult 
population and the Laplacian of the adult population and 
in the other case the situation is precisely reversed. 

In both cases there is a homogeneous nontrivial non- 
negative steady state solution, for Eqs. |(T]), ([2]) or ([3]) 
that can be written in terms of the parameters of the 
problem (3, S, /i and K, and is independent of the inter- 
action strength n: 
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(4) 



where A — \J ^ 4:Pfi — 2 5^ + fi'^. In the following sec- 
tion we will study Eqs.|[2|)and |[3]), i.e. we will analyze 
the effects of each of the new terms separately. 

We are interested in effects of temporal or seasonal 
changes in the environment. These changes will be re- 
flected in a change in the carrying capacity of the habi- 
tat. Therefore, we will consider a varying K. The ability 
of a given species to adapt to a changing environment 
is related to their mobility as well as to other quantities 
such as the birth and death rates. One way to charac- 
terize the mobility is by considering the so-called Fisher 
velocity of the species associated with traveling environ- 
ments which can naturally arise as the seasons change. 
When favorable conditions move in space and time, it is 
important to find out how the population can or cannot 
follow them, and whether there are critical velocities of 
the traveling environment which separate parameters re- 
gions in which the species survive or undergo extinction. 
To facilitate the subsequent discussion we introduce the 
concept of a refugium, a bounded domain with a high 
carrying capacity. The population of a given species can 
live within this region. Outside the refugium, the living 
conditions are too harsh for the species to survive. 



STATIC REFUGIUM 

When no temporal variations of the refugium are con- 
sidered, we find stationary profiles for both populations 
which differ from those previously found when consider- 
ing Fisher-like coupled equations. In FigsH] we display 
the profiles of both populations in three cases. All pa- 
rameters of the problem are the same in all cases, with 
the exception of k which is in (a) but, in arbitrary units, 
10 in (b) and (c). 

We can see that while the profile of the adult density 
Ma remains almost unchanged in the presence of inter- 
actions, that of the subadult density My presents inter- 
esting effects. Since the interactions tend to drive the 
subadults away from the adults, in both cases there is an 
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FIG. 1: Steady state profiles of Ma (solid)and My (dashed)in 
a bounded domain. The plots correspond to the solution to 
(a) Eqs.lH]), (b) Eqs.JJl), (c) Eqs.©. Case (a) represents no 
interaction whereas cases (b) and (c) describe the two differ- 
ent interactions considered (see text). Parameters (arbitrary 
units) are /3 = 4, ^ = 0.4, 5 = 0, k = 10. Additionally the 
environment parameter K is 10"* (practically zero) outside 
and 9 inside the refugium. 



increase of the number of individuals at the borders of the 
refugium. The effect of the interaction term is more evi- 
dent in regions where the gradient of Ma is greater. This 
fact explains the shoulder on the profile of My in Fig. 
[TJb. On the other side, the profile in Fig. [TJc presents 
the smoothest shape, consistent with a greater effective 
diffusion coefficient. 

The interaction term affects not only the shape of the 
profile of both population but also the maximum values 
attained by the population inside the refugium. We note 
that the solutions in Eq.|[4l), where n plays no role, corre- 
spond to infinite domains. When bounded domains are 
considered, Z3a, Dy^ as well as k and the domain size, 
affect the maximum value of the population. We plotted 
In Fig. [2] the behavior of the maximum attained by My, 
M™°^ , for a given size of the refugium and constant dif- 
fusion coefficients as a function of the values of k. The 
population of subadults undergoes an overall decay as k 
increases. The curves in Fig. ^ were fitted with decay- 
ing exponentials. A mono-exponential fit works in each 
case. 

An interesting effect is the possibility of maintaining 
a subadult population localized when it is bounded by a 
population of adults as if they were forming a well con- 
fining the younger mice. As the interaction terms pre- 
vent the population of subadults from moving towards 
an increasing gradient of adults, the subadult population 
remains localized until the population of adults has de- 
creased sufficiently (by diffusion) to reach a state in which 




K 

FIG. 2: Maximum value of the subadult mouse population 
My as a function of the strength k of the interaction between 
subadults and adults from Eqs. (2) (full circle, solid line ) and 
Eqs. (3) (empty circle, dashed line). The curves are the result 
of exponential fitting of the solutions obtained numerically. 



diffusion and interaction in the subadult population can 
compete and delocalize the nucleus of subadults. 



TRAVELING REFUGIUM 

Consider now a traveling refugium of constant size. 
This will allow us to test the ability of the population to 
survive while following the moving environment. We will 
find the critical velocity of the refugium Vc, above which 
survival is no longer possible. Fig. [3]displays the changes 
in the fronts as the refugium starts to move in each of 
the three cases. 

Again, while the profile of Ma (adults) remains almost 
unchanged, the profile of My (subadults) shows some ef- 
fects. Both interactions (expressed in Eqs. [2] and [3] tend 
to maintain the symmetry of the distribution of the pop- 
ulation despite the translational movement; the effect of 
Eqs. [3] is with stronger intensity as it strongly biases the 
population contrary to the movement. 

When considering only one species and the Fisher 
equation, the critical velocity is given by the Fisher ve- 
locity Vf — V Da, with D the diffusion constant and a the 
birth rate. When dealing with a set of coupled equations 
it is still possible to show that there is a critical velocity, 
associated closely with the critical velocity of the slowest 
species. 

If the new terms in Eqs. ([2]) and ^ are neglected to get 
Eqs. ([l]), within the range of the values used for the pa- 
rameter values of both populations, adult and subadult 
densities can be shown to be of the same order through- 
out the whole domain. The adult population would then 
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FIG. 3: Steady state profiles of Ma (solid)and My (dashed)in 
a traveling bounded domain with velocity « = 1.3. The plots 
correspond to the solution to (a) Eqs.([T|), (b) Eqs.((2]), (c) 
Eqs.©. With /3 = 4, = 0.4, 5 = 0, = 10 and if = 9 
inside the bubble. Parameters are in arbitrary units. 




FIG. 4: Critical velocity of the population (relative to the 
Fisher velocity) as a function of k from Eqs. (2) (full circle, 
solid line ) and Eqs. (3) (empty circle, dashed line). The 
curves are exponential fitting of the data. 

be described by 

^ = D,V2M, + AiM,-<5M,-9^. (5) 
at K 

This would mean that we can associate a Fisher veloc- 
ity, Vf — ^ Dain — S) with the population Ma. While it 
appears that a similar argument might be used for My, 
this is not true because the idea that both populations 
are of the same order is no longer valid in that case when 
the velocity of the refugium is above the Fisher veloc- 



FIG. 5: Limiting maximum values M^"^^ (squares) and 
M™"^" (circles). Full symbols correspond to Eqs. (2) while 
empty symbols correspond to Eqs. (3). Curves are exponen- 
tial fittings of the data 

ity for Ma. The adults are slower than subadults. It 
is therefore that the lower Fisher velocity appropriately 
describes the critical situation. This explains the numer- 
ical results that show that when considering Eqs|(l]), the 
critical velocity of the entire population Vc is almost the 
same as the critical velocity that can be obtained from 
Eq.©. 

We display in Fig. |4] interaction effects on the critical 
velocity Vc by plotting Vc/vf (where Vf is the Fisher ve- 
locity of the adult mice in the absence of interactions) as 
a function of the interaction strength k for the two kinds 
of interaction. We find that the numerical solutions may 
be fitted excellently by saturating exponentials. We ob- 
serve that in both cases there is an increase in Vc as n 
increases. The effect is much more evident for EqlJ] as 
was expected from the shape presented by its profile in 
Fig. S). 



BREATHING REFUGIUM 

Seasonal variations of the relevant parameters may cor- 
respond to the refugium being centered in a static point 
with its size changing in an oscillating way. We call this 
case a breathing refugium. To show the new features due 
to the effect of the interaction terms, we plot the maxi- 
mum of the subadult density My as a function of time. If 
there are no interactions between subadults and adults, 
the maximum suffers negligible periodic variations. On 
the other hand, in the presence of interactions, the max- 
imum oscillates between two values displayed in Fig. [5] 
by squares and circles, respectively. We call these two 
limiting maximum values M™"^"'" and M""*^~. The am- 
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FIG. 6: Relative amplitude of the temporal oscillations of the 
maximum value of My{x,t) as a function of k from Eqs. (2 
)(full circle, solid line ) and Eqs. (3) (empty circle, dashed 
line). The curves are exponential fitting of the data. 
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FIG. 7: Comparative intensity of each effect on the border of 
the moving front of subadult population. The effect associ- 
ated to Eqs. (I2|) is in solid line, while the one corresponding 
to Eqs. ^ is in dashed line. The units are arbitrary 



plitude of this variation increases in a very apparent way 
as K grows. This is shown in Fig. [6l where we plot the 

relative difference a — " ,,m.a^l , which is the ratio 

of the difference between the limiting maximum values 
to the greater of these maximum values. 



CONCLUSIONS 

We have shown the emergence of some interesting ef- 
fects due to the inclusion of interaction among differ- 
ent groups in a population. We have also compared our 
new results with those already known from simple Fisher 
equations with no interactions. Among the new findings 
are the change in the shape of the steady and travel- 
ing profiles (Figs. [T] and [H, the increase of the criti- 
cal velocity of the populations (Fig. |4]) with interaction 
strength, and the oscillatory behavior of the population 
profiles when seasonal changes in the environment are 
considered (Figs. [5]and[6l). Some of these results can be 
understood easily in the context of Fig. [H When a trav- 
ehng refugium is considered, the steady fronts displayed 
in Fig. [T] change. The traveling profile are asymmetric, 
showing a depletion of the concentration of individuals 
at the head (right) of the front. In the presence of the 
interaction, the subadult population tends to move to- 
wards areas less populated with adults. This explains 
the changes in the traveling profiles presented by both 
population in Figs. [SK and [Sb. The effect of the term 
in Eqs. [2] is stronger than that of the one considered in 
EqslSl Indeed, we observe that not only the population 
of subadults moves to the right in a more evident way but 
the change in the critical velocity is also more intense. 

The shape presented by the profiles associated to each 
of the equations can be understood by qualitative ar- 
guments. As mentioned before, a comparison of the two 
new interaction expressions, shows that they share a com- 
mon term proportional to the gradient of each of the 
classes V(Ma)(VMy). The other terms, Mj^V^Mq and 
MaV^Mjy respectively, promote a fiux of the subadult 
population towards less populated regions, i. e. to the 
right of the traveling front. This explains the shift of 
the profile of subadults. The difference in the resulting 
profile can be understood when analyzing the difference 
in intensity of ech term, as shown in Fig. [71 where we 
compare the intensity of each of the effect when consider- 
ing a stationary profile. The oscillations found when the 
refugium is breathing are due to the fact that a system 
[h]iescribed by any of the Eqs. |[2l) or ([3]) is much more 
sensitive to changes in the size of the refugium, as can be 
observed by the values displayed in Fig. [21 where with- 
out loss of generality, only one size of the refugium was 
analyzed. 

This work was supported in part by the NSF under 
grant no. INT-0336343 and by NSF/NIH Ecology of In- 
fectious Diseases under grant no. EF-0326757. 
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